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ABSTRACT 

Plasma  jet  generators  have  been  designed  and  tested  which  used  an  explosive  driver  and 
shocktube  with  a  rectangular  cross  section  that  optimize  the  flow  velocity  and  electrical 
conductivity.  The  latest  in  a  series  of  designs  has  been  tested  using  a  reactive  load  to  diagnose 
the  electrical  properties  of  the  MHD  generator/electromagnet  combination.  The  results  of 
these  tests  indicate  that  the  plasma  jet/MHD  generator  design  does  generate  a  flow  velocity 
greater  than  25  km/s  and  produces  several  gigawatts  of  pulsed  power  in  a  very  small  package 
size.  A  larger,  new  generator  design  is  also  presented. 


INTRODUCTION 


There  have  been  many  studies  and  experiments  conducted  over  the  past  several  decades  (Refs.  1,  2,  and  3)  on 
explosively-driven  MHD  generators  (EDMHDG).  A  recent  study  of  such  devices  (Ref.  4)  demonstrated  an 
explosive  driver  design  that  produced  1.8  GW  of  pulsed  power  using  argon  gas  at  a  velocity  of  25  km/s  with  a 
conductivity  of  18,500  S/m.  These  early  experiments  were  done  using  resistive  loads.  In  the  current  effort,  the 
latest  available  design  was  investigated  for  its  electrical  properties  using  a  reactive  (LC)  load.  Finally,  a  larger 
scale  design  was  produced  to  reach  higher  voltages  and  output  power  levels. 


MHD  GENERATOR  DESIGN 


Figure  1  depicts  the  design  of  the  EDMHDG  used  in  the  experiments  reported  here.  Two  pounds  of  explosive 
(COMP  C4)  are  detonated  from  the  left  hand  side  using  a  exploding  bridgewire  (EBW)  Reynolds  Electronics  RP-1 
detonator.  The  explosive  compresses  via  a  shaped  charge  liner,  heats,  and  ionizes  the  argon  gas  in  the  adjoining 
cavity  initially  at  200  psi  pressure.  The  mylar  diaphragm  bursts  with  the  force  of  the  shock  wave  and  allows  the 
ionized  argon  to  flow  into  the  evacuated  lexan  shock  tube,  or  MHD  channel.  An  electromagnet  is  energized  by  the 
discharge  of  a  capacitor  bank  so  as  to  produce  peak  field  during  the  main  flow  of  Argon  plasma  in  the  channel. 
The  subsequent  U  x  B  force  (where  U  is  the  flow  velocity  of  the  plasma  and  B  is  the  external  magnetic  field)  on 
the  free  electrons  in  the  plasma  causes  charge  to  be  collected  at  the  terminals  to  which  the  generator  load  is 
connected. 


*  Sandia  Contractor  to  this  work  was  supported  by  the  US  Department  of  Energy  under  DE-AC04-94AL85000. 
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MHD  GENERATOR  DESIGN  EXPECTATIONS 


The  open  circuit  voltage  developed  at  the  electrodes  in  the  MHD  channel  is  given  by: 

VM  =  UBd  (1) 

where,  U  =  flow  speed  in  m/s,  B  =  magnetic  field  strength  in  Tesla,  and  d  =  distance  between  the  pickup 
electrodes  in  meters.  The  power  delivered  to  the  load  for  the  MHD  generator  for  the  low  magnetic  Reynolds 
number  case  can  be  estimated  for  a  rectangular  channel  from  (Refs.  1,4): 

P  =  K(l-K)a-(UB)2(Vol),  for  Rm*  1  (2a) 

„  2B2  Udw 

P  =  - ,  for  Rm»  1  (2b) 

JUo 

where 

P  =  power  in  Watts 

K  =  generator  load  parameter  (dimensionless) 

ct  =  electrical  conductivity  of  the  plasma  in  Siemens/m 
p0  =  permeability  of  free  space  =  4n  x  10'7  H/m 

w  =  inside  width  of  rectangular  MHD  channel  in  meters 
Vol  =  Volume  of  plasma  slug  in  rectangular  channel  in  m3 

The  generator  load  parameter  is  defined  as  the  ratio  of  the  generator  output  voltage  under  load  to  the  open  circuit 
generator  voltage.  Maximum  power  output  is  delivered  to  the  load  when  K  =  0.5.  The  magnetic  Reynolds 
number,  Rm  =  p0aUd,  is  a  measure  of  the  magnetic  field  trapped  in  the  flow  of  the  plasma.  When  Rm  -  one  can 
assume  no  magnetic  field  is  trapped  near  the  boundaries  of  the  channel.  Note  that  for  the  very  high  magnetic 
Reynolds  number  case,  the  power  delivered  by  the  generator  is  independent  of  the  conductivity  of  the  plasma  and 
the  length  of  the  plasma  slug.  Also,  the  power  only  depends  linearly  on  the  flow  velocity.  As  will  be  seen  below, 
the  MHD  generators  in  these  experiments  have  Rm  greater  than  10.  For  this  type  of  generator,  the  Hall  effect  is 
important,  the  current  in  the  plasma  reaches  a  limit  of  2B0w/p0 ,  and  the  most  important  control  parameters  are  the 
flow  velocity  and  external  magnetic  field. 

The  EDMHD  Generator  design  #6  (Figure  1)  configuration  was  modeled  using  the  CTH  hydrocode  (Ref.  5) 
developed  at  Sandia  National  Laboratories.  The  CTH  code  predicted  flow  parameters  at  the  center  of  the  magnetic 
field  produced  by  the  electromagnet  (9  inches  from  the  shocktube  entrance  and  diaphragm).  The  key  results  of 
CTH  code  calculations  for  design  #6  using  2  lbs  of  explosive,  initial  argon  gas  pressure  of  200  psi,  and  a  shocktube 
air  pressure  of  0.01  psi  were:  argon  density  =  44.5  kg/m3,  static  pressure  behind  shock  =  2.7  kbars,  plasma 
temperature  =  23,000  K,  plasma  flow  velocity  =  25.5  km/s,  and  plasma  flow  arrival  time  =  57.7  ps. 

In  order  to  estimate  the  power,  the  electrical  conductivity  must  be  determined.  For  the  case  of  a  singly  ionized  gas, 
the  Spitzer  conductivity  can  be  approximated  by1 : 


0.943 ^ 

l7Wi 

8  me 

e2ln 

T  3(kT)% 

1  See  paper  in  Shuler  (Ref.  8) ,  "Ionization  Behind  Shock  Waves,"  by  J.  D.  Teare,  p.  222. 
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In  this  expression,  k  =  Boltzmann  constant  =  1.38  x  10'16  erg/K,  T  =  electron  temperature  in  Kelvin,  me  =  mass  of 
electron  =  9. 1 1  x  10'28  g,  and  e  =  electronic  charge  =  4.8  x  10'10  esu. 

The  pressure  behind  the  shock  in  the  argon  gas  was  calculated  to  be  2.7  kbars  or  2.7  x  108  Pa,  gas  density  of  44.5 
kg/m3>  and  the  flow  velocity  is  25.5  km/s.  The  conductivity  can  now  be  estimated.  The  correct  temperature  to  use 
in  eq.  (3)  is  the  electron  temperature.  With  the  argon  gas  speed  (atom  velocity)  estimated  at  25.5  km/s,  a 
reasonable  value  for  electron  temperature2  is  14,000  K  if  a  two  percent  degree  of  ionization  is  assumed.  Under  the 
above  conditions,  the  electron  density,  ne,  is  1.34  x  1019  cm'3.  Thus,  the  electrical  conductivity  of  the  argon  gas 
flow  is  estimated  to  be  1.66  x  1014  s'1,  or  18,400  Siemens/m. 

The  experiment  discussed  here  has  been  designed  to  include  a  magnetic  field  in  the  plasma  flow  channel,  or 
shocktube,  of  about  9.4  T  using  two  coils  (one  on  each  side  of  the  shocktube)  consisting  of  13  turns  of  litz  wire 
cable  (Ref.  6).  The  energy  for  these  coils  is  supplied  by  a  519  pF  capacitor  bank  charged  to  10.43  kV  for  28.2  kJ  of 
stored  energy. 

The  expected  results  of  this  experiment  can  now  be  determined  for  the  0.5  x  2  inch  cross  section  shocktube 
(electrodes  are  positioned  across  the  2  inch  gap,  thus  d=  0.0508  m).  From  eq.  (1),  the  open  circuit  voltage 
generated  by  the  EDMHDG  is  estimated  to  be  12.2  kV.  The  magnetic  Reynolds  number  is  30.  Therefore,  from  eq. 
(2b),  the  maximum  power  delivered  to  the  load  is  2.3  GW  .  This,  of  course,  implies  a  peak  current,  Im,  through  the 
generator  of  189  kA. 


TEST  SETUP  AND  PROCEDURES 


The  MHD  generator  depicted  in  Fig.  1  (design  #6)  was  tested  using  a  capacitive  load  on  two  separate  tests.  The 

load  was  designed  to  match  the  impedance  of 
the  generator.  An  equivalent  circuit  model 
along  with  laboratory  pulser  testing  of  the 
load  were  used  to  determine  the  sensitivity  of 
the  experiment  to  the  dynamic  impedance  of 
the  generator.  The  main  objective  of  these 
tests  was  to  determine  a  more  exact 
measurement  of  the  dynamic  impedance  of 
design  #6.  The  load  circuit  was  instrumented 
with  voltage  dividers  and  Rogowski  coils.  In 
addition,  the  MHD  generator  was 
instrumented  with  sensors  to  measure  the 
open  circuit  voltage,  time  of  plasma  arrival 
along  the  shocktube,  electromagnetic  coil 
performance,  eddy  current,  etc.  Table  1 
provides  a  list  of  all  these  diagnostics  and  Ref. 
7  explains  their  use  and  operation  in  more 
detail. 


Figure  1.  EDMHD  Generator  Design  #6 


The  procedures  followed  during  the  two  tests  were  as  follows:  the  capacitor  bank  was  charged  to  a  nominal  10 
kV.  Once  full  charge  was  reached,  a  signal  from  the  fire  set  triggered  the  bank  to  discharge.  Then,  a  delay 
generator  caused  the  detonator  to  fire  94  ps  after  the  fire  set  signal.  This  delay  was  timed  to  allow  for  maximum 
plasma  flow  in  the  channel  at  peak  magnetic  field.  Data  was  recorded  from  the  probes  for  about  7  ps  pulse 
duration. 


2In  Shuler,  "Ionization  Behind  Shock  Waves,"  by  J.  D.  Teare,  p.  229  and  the  plot  on  p.  267. 
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EXPERIMENT  RESULTS 


A  circuit  simulation  using  MICROCAP  was  conducted  to  determine  the  variation  of  load  voltage  and  current  as  a 
function  of  generator  total,  dynamic  impedance.  The  model  constructed  to  simulate  the  tests  is  shown  in  fig.  2. 
The  results  of  the  circuit  simulation  are  depicted  in  figs.  3  and  4. 

Table  I.  EDMHD  Generator  Probes  and  Related  Parameters 


No. 

LD./Data 

Channel 

Type 

E 

Related 

Parameter 

Comment 

1 

OCVP#l 

Voltage 

7.85 

Vn„R.a 

Open  circuit  voltage  -  high 

2 

OCVP  #1A 

Voltage 

7.85 

V„„,  R,  a 

Open  circuit  voltage  -  low 

3 

PIN  SW  #5 

Arrival  time 

30.75 

U 

Ionization  pin/plasma  arrival 

4 

BW  current 

Current 

N/A 

N/A 

Detonator  bridgewire  current 

5 

Current 

20.75 

a 

Eddy  coil  current  -  high 

6 

ttr.mfsmm 

Current 

20.75 

cr 

Eddy  coil  current  -  low 

7 

Current/time 

20.75 

I  dot 

Differential  current  -  high 

8 

Eddy  di/dt  #1A 

Current/time 

20.75 

I  dot 

Differential  current  -  low 

9 

Volt.  Div.  #1A 

9.75 

V 

Load  voltage  -  high 

10 

Volt.  Div.  #2A 

9.75 

V 

Load  voltage  -  low 

11 

ISOE9H 

Current 

9.75 

I 

Load  current  -  high 

12 

Current 

9.75 

I 

Load  current  -  low 

13 

Volt.  Div.  #1 

9.75 

V 

Load  voltage  -  high 

14 

Volt.  Div.  #2 

9.75 

V 

Load  voltage  -  low 

15 

Rog.  #1 

Current 

9.75 

I 

Load  current  -  high 

16 

Rog.  #2 

Current 

9.75 

I 

Load  current  -  low 

17 

PIN  SW  #1 

Arrival  time 

1.75 

u 

Ionization  pin/plasma  arrival 

18 

PIN  SW  #2 

Arrival  time 

4.00 

u 

Ionization  pin/plasma  arrival 

19 

PIN  SW  #3 

Arrival  time 

22.75 

u 

Ionization  pin/plasma  arrival 

20 

PIN  SW  #4 

Arrival  time 

26.75 

u 

Ionization  pin/plasma  arrival 

21 

EM  Coil  Rog. 

Current 

N/A 

B 

EM  coil  current 

22 

EM  Coil  Volt. 
Divider 

Voltage 

N/A 

B 

EM  coil  voltage 

^Measured  from  shocktube  entrance 


The  simulation  results  indicate  that  the  designed  load  provides  ample  sensitivity  to  measure  the  impedance  of  the 
generator  from  10  mohms  to  over  100  mohms.  Figure  5  contains  the  load  voltage  data  for  each  experiment. 

Figure  6  shows  the  data  from  one  of  the  open  circuit  voltage  probes  for  both  tests.  As  can  be  seen  from  the  data, 
the  peak  voltages  measured  were  considerably  higher  than  the  simple  estimates  discussed  earlier.  For  example,  the 
open  circuit  voltage  for  both  shots  exceeded  20  kV  and  reached  over  30  kV  for  Test  #1.  Assuming  a  simple 
matched  load,  the  expected  voltages  on  the  load  would  be  one  half  the  open  circuit  voltage.  Figure  5  shows  a  peak 
voltage  that  does  approach  15  kV  but  does  not  quite  reach  it.  In  analyzing  the  waveform  for  the  load  data  on  the 
first  test,  the  best  match  that  could  be  made  to  the  simulation  curves  indicates  a  generator  impedance  of  130 
mohms.  In  addition,  the  PIN  switch  data  for  both  tests  suggests  a  particle  velocity  of  20  km/s  which  is  20  percent 
lower  than  the  CTH  calculated  value.  The  peak  power  delivered  was  1.8  GW,  22  percent  less  than  predicted. 

It  should  be  noted  that  electrical  breakdown  probably  occurred  on  both  tests  before  a  full  pulse  of  energy  was 
delivered.  The  breakdown  is  clearly  evident  on  the  second  test  and  occurred  before  the  peak  voltage  was  reached. 
The  data  from  load  on  the  first  test  also  shows  evidence  of  electrical  breakdown  at,  or  just  following,  the  peak 
signal.  It  is  believed  that  a  voltage  probe  cable  breakdown  occurred  in  the  second  instance.  The  origin  of  the 
breakdown  for  the  first  test  has  not  been  definitely  established. 
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Figure  2.  MICROCAP  Circuit  Simulation 


Figure  3.  Circuit  Simulation  Results  -  Current  vs.  time 
with  Load  Impedance  as  Free  Parameter 


Load  Vol  t  age 


Figure  4.  Circuit  Simulation  Results  -  Voltage  vs.  time  with  Figure  5.  Load  Voltage  Test  Results 

Load  Impedance  as  Free  Parameter 


Open  Circuit  Voltage  Probe 


Figure  6.  Open  Circuit  Voltage  Test  Results 


Figure  7.  EDMHD  Generator  Design  #7 
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LARGER.  SCALED  EDMHD  GENERATOR  DESIGN 


A  new  design  has  been  completed  and  hardware  has  been  built  and  partially  tested.  The  purpose  of  the  new  design 
is  to  determine  if  a  much  larger  rectangular  shocktube  will  allow  higher  voltages  with  similar  currents  (thus,  at 
higher  impedance)  to  be  delivered  to  a  capacitive  load.  To  test  this  concept,  an  EDMHD  Generator  with  a 
shocktube  that  is  12  in.  x  0.5  in.  x  36  in.  long  has  been  constructed.  Figure  7  shows  the  basic  design  of  this  larger 
generator.  A  new  electromagnet  design,  depicted  in  the  figure  as  the  U-shaped  object,  was  required  to  provide  a 
field  of  12  - 15  T  within  the  large  volume  of  the  shocktube.  This  coil,  each  half  made  of  4130  steel  plates  welded 
together  to  make  10  turns  and  connected  so  as  to  have  additive  fields  between  each  half,  has  been  tested  up  to 
energies  of  1.5  MJ.  The  peak  field  along  the  axis  of  the  channel  was  measured  at  12.2  T  for  this  case.  The  volume 
that  this  field  is  maintained  at  its  peak  is  contained  within  the  3  inch  by  13  inch  oval  slot  with  a  spacing  of  2.25 
inches. 

The  following  performance  estimates  are  determined  from  eqs.  (1)  -  (3)  and  the  experience  and  modeling 
performed  on  the  previous  design  and  using:  U=  25.5  km/s,  electron  temperature  =  14,000  K,  ct=  18,400 
Siemens/m,  Then,  Voc  =  95  kV,  Rm  =  180,  P=  23  GW,  andlm=  242  kA.  Note  that  by  scaling  the  shocktube 
by  a  factor  of  six  in  electrode  separation  increases  the  open  circuit  voltage  by  the  same  factor.  Since  the 
conductivities  are  about  the  same  for  both  cases  with  a  high  magnetic  Reynolds  number,  the  power  output  is  driven 
primarily  by  the  size  of  the  magnetic  field.  Since  the  power  increases  as  the  square  of  the  magnetic  field,  a  small 
increase  in  electromagnetic  coil  current  will  produce  a  significant  amount  of  more  power.  This  increase  will  have 
to  be  weighed  against  the  problems  associated  with  producing  such  high  fields  in  a  larger  volume  as  the  generator 
increases  in  size.  This  latest  design  of  the  EDMHD  Generator  has  been  built  and  will  be  tested  experimentally  to 
determine  the  effects  of  the  many  nonlinear  properties  associated  with  the  MHD  flow. 


REFERENCES 


1.  Jones,  Jr.,  M.  S.  and  Blackman,  V.  H.,  "Parametric  Studies  of  Explosive  Driven  MHD  Power  Generators," 
MHD  Electrical  Generation,  Vol.  2,  1964,  p.  803. 

2.  Jones,  Jr.,  M.  S.,  and  McKinnon,  Jr.,  C.  N.,  "Explosive  Driven  Linear  MHD  Generators,"  Conference  on 
Megagauss  Magnetic  Field  Generation  by  Explosives  and  Related  Experiments,  Knoepfel,  H.  and  Herlack,  F. 
(eds.),  European  Atomic  Energy  Community:  Brussels,  1966,  p.  349. 

3.  Sharma,  S.  P.,  Schmidt,  H.  J.,  et.  al.  "Feasibility  of  an  Optically  Pumped  Iodine  Excimer  .342  mm  Laser  Using 
Explosive  MHD  Generator  Plasma,"  Megagauss  Technology  and  Pulsed  Power  Applications,  Fowler,  C.  M., 
Caird,  R.  S.,  and  Erickson,  D.  J.  (eds.),  Plenum  Press:  New  York,  1987,  p.  757. 

4.  Vigil,  M.  G.,  Duggins,  B.  D.,  Explosively-Driven  Magnetohydrodynamic  Generator:  Phase  I,  Sandia 
National  Laboratories,  SAND94-0045,  Aug.  1994. 

5.  McGlaun,  J.  M.,  Thompson,  et.  al.,  A  Brief  Description  of  the  Three-Dimensional  Shock  Wave  Physics 
Code  CTH,  Sandia  National  Laboratories,  SAND89-0607,  1990. 

6.  Kaye,  R.  J.,  Cnare,  E.  C.,  et.  al.,  "Design  and  Performance  of  Sandia's  Contactless  Coilgun  for  50  mm 
Projectiles,"  IEEE  Trans,  on  Magnetics,  29(1),  pp.  691-695,  Jan,  1993. 

7.  Vigil,  M.  G.,  Duggins,  B.  D.,  Diagnostics  Instrumentation  Systems  Development:  Explosively-Driven 
Magnetohydrodynamic  Generator  Project;  Phase  I,  Sandia  National  Laboratories,  SAND94-0046,  February, 
1994. 

8.  Shuler,  K.  E.  and  Fenn,  J.  B.  (eds.),  Ionization  in  High-Temperature  Gases,  Academic  Press:  New  York, 
1963. 


1073 


